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Using density functional theory, we examine a nitrogen-doped anionic golden cage (NAu^~ 6 ). For 
the exohedral adsorption that is more stable than the endohedral doping, the bridge and hollow sites 
have larger binding energies than the atop sites by ^1 eV. When the N atom is adsorbed on the 
cage, electrons are transferred from AuJ" 6 to the N atom. The transition between the exohedral and 
endohedral adsorption may occur thermally through a bridge site. In the infrared active vibrational 
spectra, exohedral doping causes greater intensities at higher frequencies than endohedral doping. 



PACS numbers: 61.48.C, 68.43.Bc, 73.22.-f 

INTRODUCTION 

Since their discovery, carbon fullerene molecules such 
as Ceo, C70, and Cgo have attracted considerable at- 
tention because of their remarkable electronic, optical, 
and electrochemical properties QUI]- Interestingly, it 
was proved experimentally that carbon fullerenes have 
cousins made of metal, in particular, golden fullerenes. 
In 2006, S. Bulusu and colleagues showed experimental 
and theoretical evidence of anionic hollow golden cages, 
Au~ (n = 16 — 18) [5]. They compared experimental pho- 
toelectron spectra of the anionic structures with those 
simulated by density functional theory (DFT) [6]. 

Au^ 6 is the smallest golden cage with an empty inte- 
rior. It has slightly distorted tetrahedral symmetry (T^), 
with an inner diameter of ^5.5 A. Thus, it is possible to 
endohedrally dope the golden cage in analogy to the car- 
bon fullerenes. Numerous theoretical and experimental 
studies of golden cages adsorbed with metal atoms have 
been reported [THTO]. Similar to carbon fullerene deriva- 
tives, this kind of golden cage can be a host to form var- 
ious derivatives by adsorbing various atoms, molecules 
and radicals. Recently, due to their high reactivity, gold 
nanoclusters have been explored for their potential use 
for nanocatalysts [TTHT4] . 

In this paper, we report our study on a golden cage 
anion, Auj~ 6 , doped with a nitrogen atom which is a non- 
metal atom [15]. Since oxidation is an important issue 
for gold nanostructures because of their high reactivity, 
we investigate the binding characteristics of the adsor- 
bate atom on Auj~ 6 . Nitrogen is an element with a large 
electron affinity, and so the adsorbed N atom may ac- 
cept electrons from the golden fullerene. Going a step 
further, it may help to understand adsorption properties 
of molecules containing nitrogen such as N2, NO, NO2, 



and NH3, onto the golden buckyball. We also investigate 
the transition between the most favorable exohedral and 
endohedral adsorption configurations. Finally, we show 
the infrared (IR) active vibrational spectra of our mod- 
els, which may help identify N-adsorbed Au^ 6 in future 
experiments. 



COMPUTATIONAL DETAILS 

To study the binding properties of a nitrogen atom on 
Auj~ 6 , we performed first-principles calculations within 
DFT with spin polarization implemented in the DMol 3 
package [16]. We used the generalized gradient approxi- 
mation with the Perdew-Burke-Ernzerhof functional [17] 
to describe the exchange-correlation functional. To ex- 
pand the electronic wavefunctions, a double- numerical 
polarized basis set was chosen with a real space cutoff of 
6.5 A. The scalar relativistic effects were included in the 
all-electron calculation. For more precise computations, 
we chose an octupole scheme for the multipolar fitting 
procedure on the charge density and a fine grid scheme 
for numerical integration [T8] . 

The cluster geometry was optimized by the Broyden- 
Fletcher-Goldfarb-Shanno algorithm [T9H23] without 
symmetry constraints until the total energy change was 
converged to 10 -5 Ha in the self-consistent loop and the 
atomic forces were smaller than 4 x 10 -4 Ha A -1 . The 
optimized geometries of NAuj~ 6 were determined from 
the total energy minimum, and were further verified by 
the absence of imaginary frequencies in the harmonic 
frequency calculations. Using a nudged elastic band 
(NEB) method [24], we searched for the minimum-energy 
path in geometrical phase space connecting two stable 
structures. Including the initial and final configurations, 
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FIG. 1: Atomic model structures of the nitrogen atom adsorbed on Au^~ 6 . Ten adsorption sites of the adsorbate are considered: 
three different triangular sites, TR1, TR2, and TR3; three dissimilar bridge sites, BR1, BR2, and BR3; two inequivalent atop 
sites, ATI and AT2; and two endohedral doping sites, CTR and In-BR3. For the detailed classification of these adsorption 
sites, see the text. 



twenty replicas were chosen to construct the NEB. For 
a trajectory of the minimum-energy path between reac- 
tants and products during the atomic encapsulation, lin- 
ear synchronous transit and quadratic synchronous tran- 
sit (LST/QST) calculations [25, 26 were also performed 
with conjugate gradient minimization [27 . Finally, IR- 
active vibrational spectra were evaluated for our models 
by calculating the dynamical matrix [28] and the Born 
effective charge tensor [29] . 

The golden cage Au^~ 6 is a truncated tetrahedron with 
tetrahedral symmetry. As in our previous study [15], we 
considered various endohedral and exohedral adsorption 
sites for the initial configurations, as shown in Figure [T] 
A truncated tetrahedron is composed of four truncated 
triangular facets (defined as TR1) [15] and four hexag- 
onal facets. In particular, each hexagonal facet in Auj~ 6 
consists of six triangles as shown in Figure [T] They are 
geometrically different from TR1, and are classified into 
two types, TR2 and TR3. TR2 shares an edge with 
an adjacent TR1 and two edges with two TR3 triangles 
in the same hexagonal facet, whereas TR3 has a com- 
mon junction with another TR3 located on a neighboring 
hexagonal facet, and two junctions with two TR2 trian- 
gles in the same hexagonal facet. Three TR2 and three 
TR3 triangles are arranged alternately on each hexag- 
onal facet. Edges shared by two neighboring triangles 
or bridges between two Au atoms can be classified into 
three types, BR1, BR2, and BR3. BR1 is shared by TR1 
and TR2; BR2 is done by TR2 and TR3 in the same 
hexagonal facet; and BR3 by two TR3 triangles, one in a 
hexagonal facet and the other in a neighboring hexagonal 
facet. Besides, two types of vertices are considered (atop 
sites): ATI and AT2. ATI represents twelve equivalent 



vertices forming four TR1 triangles, and the remaining 
four vertices located at the center of each hexagonal facet 
fall into the other type, AT2. Overall, there exist eight 
different adsorption sites on the outer surface of Au^ 6 as 
shown in Figure [I] 

We also considered the encapsulation inside Au^ 6 
as "adsorption" sites for endohedral doping of the N 
atom [15 . Together with the center (CTR) of the Auj~ 6 
cluster, eight more off-centered encapsulation configura- 
tions named as "In-XXn" corresponding to the eight ex- 
ohedral sites (XXn, XX=TR, BR, or AT). Our DFT cal- 
culations revealed that only CTR and In-BR3 sites are 
stable among all endohedral doping sites, and the other 
seven adsorption configurations are unstable and tend to 
change into either In-BR3 or CTR during structural re- 
laxation. 

RESULTS AND DISCUSSION 

For each configuration shown in Figure [I] we obtained 
the equilibrium geometry of NAu^~ 6 and its binding en- 
ergy (E\y) defined by 

E b = £[Aur 6 ] + E[N] - £[NAur 6 ], 

where £?[NAuj~ 6 ] and i?[Auj~ 6 ] are the total energies of the 
golden cage with and without an N atom, respectively, 
and E[N] represents the energy of an isolated N atom. All 
the calculated energy values are summarized in Table [I] 
Here, a positive (negative) sign in means an exother- 
mic (endothermic) process. We found that two trian- 
gular sites (TR1 and TR2) and one bridge site (BR3) 
are favored by the N atom. TR1 is the most preferred 



adsorption site with a binding energy of 2.95 eV, and 
BR3 (TR2) has a binding energy of 2.90 eV (2.80 eV). 
The binding energy for the TR3 site is not as large as 
those for TR1 and TR2, as listed in Table |TJ By con- 
trast, the atop sites (ATI and AT2) are less preferable 
(«1.7 eV at ATI and wl.l eV at AT2). The N atom 
that was initially located at BR1 spontaneously moved 
to TR2 during structural relaxation, implying that BR1 
is an unstable adsorption site. 

We also calculated the encapsulation energy, which is 
regarded as an energy gain due to endohedral doping. 
For the CTR site, the encapsulation energy for the ni- 
trogen atom is 0.58 eV, which is much smaller than the 
binding energies on all the exohedral adsorption sites. In- 
triguingly, the only stable off-centered adsorption site, In- 
BR3, shows much bigger encapsulation energy (1.95 eV) 
than CTR. As shown in Figure [T] in the In-BR3 case, 
the endoherally-doped N atom appears to elongate the 
Au-Au bond corresponding to BR3 and tend to pull in 
two Au atoms corresponding to AT2, compared to its 
exohedral counterpart, BR3. 

Using the charge density difference, we investigated 
the charge redistribution induced by the N atom ad- 
sorption. Figures [2^i-c displays the charge density dif- 
ferences of TR1, the most stable exohedral adsorption 
site, and CTR and In-BR3, two stable endohedral ad- 
sorption sites. The charge differences reveal dumbbell- 
like shapes with two lobes at the positions of adjacent 
gold atoms. The charge transfer, AQ between the adsor- 
bate and the golden fullerene was also calculated using 
the Mulliken population analysis. As listed in Table [T[ 
electrons are donated from Auj~ 6 to the N atom for all 
the configurations. CTR shows the smallest amount of 
charge transfer, whereas In-BR3 does the largest charge 
transfer. We also analyzed the local magnetic moments 
of our model systems. It was found that among the ad- 
sorption configurations we considered, only CTR has the 
spin magnetic moment of ^2.1 where \±b is the Bohr 
magneton 5.788 x 10" 5 eV T" 1 , where T denotes 



TABLE I: The binding energy (E h ), the HOMO-LUMO gap 
(£?gap), the distance between the nitrogen atom and its near- 
est gold atoms in the golden cage (cZau-n), and the amount of 
electronic charge transferred from Au^ 6 to the adsorbate atom 
(AQ) at the different adsorption sites shown in Figure IT] 



N on Au 16 


configuration 


E h (eV) 


E gap (eV) d A 


a -N (A) AQ (e) 


TR1 


2.95 


0.14 


1.98 


-0.513 


TR2 


2.80 


0.83 


2.03 


-0.477 


TR3 


1.73 


0.13 


2.36 


-0.486 


BR2 


2.57 


0.32 


1.94 


-0.404 


BR3 


2.90 


0.64 


1.91 


-0.439 


ATI 


1.68 


0.57 


1.84 


-0.353 


AT2 


1.12 


0.68 


1.87 


-0.323 


CTR 


0.58 


0.69 


2.20 


-0.142 


In-BR3 


1.95 


0.26 


2.05 


-0.630 



3 




-0.51 3e -0.1 42e 




-0.630e 2.106|j B 

FIG. 2: The isosurface plots of charge density difference for 
(a) the TR1, (b) the CTR, and (c) the In-BR3 structures. The 
electron accumulation is represented in red, and the depletion 
in blue. The negative value shown at the bottom of each 
structure represents the amount of donated electrons from 
the golden cage to the N atom, which evaluated using the 
Mulliken population analysis, (d) The spin density plot of 
the CTR structure. Its net magnetic moment is calculated to 
be 2.106 /is, where \±b is the Bohr magneton. The isovalue 
is ±0.005 e A -3 . 



the tesla, the unit of magnetic flux density). Figure 
shows the spin density plot of CTR, where the majority 
spin density is distributed only at Au atoms correspond- 
ing to AT2 as well as at the N atom at the center. The 
minority spin density does not appear with the isovalue 
of ±0.005 e A" 3 . 

Next, we move on to the energy eigenvalues and the 
corresponding molecular orbitals (eigenfunctions). Fig- 
ure [3^i displays the highest occupied molecular orbitals 
(HOMOs) and lowest unoccupied molecular orbitals (LU- 
MOs) of the three configurations shown in Figures [2^l— c 
as well as those of the bare golden cage for comparison. 
Obviously, we can find the antibonding character of N 2p 
orbitals with Au 5d orbitals. The adsorption of the nitro- 
gen atom causes the breaking of the tetrahedral symme- 
try, and some degeneracy is lifted resulting in level split- 
ting as shown in Figure ^jp. The HOMO-LUMO energy 
gaps are calculated to be 0.14 eV, 0.69 eV, and 0.26 eV 
for TR1, CTR, and In-BR3, respectively, all of which 
are much smaller than that (1.37 eV) of the bare Auj~ 6 . 
Because of its much smaller HOMO-LUMO gap, the ex- 
ohedrally doped complex TR1 may be chemically more 
reactive than the endohedrally doped configurations. 
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FIG. 3: (a) HOMO and LUMO of the bare Au~ cage, the 
most stable exohedrally adsorbed structure (TR1) of NAuJ~ 6 , 
and the N-encapsulated structures (CTR and In-BR3) of 
N@Au5" 6 . The isovalue is ±0.02 e A -3 , (b) Energy levels near 
the HOMOs and the LUMOs of the four structures given in 
(a). The numbers written next to the energy levels represent 
the degeneracy of respective levels. The HOMO level of every 
structure is set to zero. For each configuration, the HOMO- 
LUMO energy gap is given at the bottom of the graph. 



To evaluate the activation energy barrier associated 
with nitrogen atom encapsulation, which plays a piv- 
otal role in the reaction rate of the dopant insertion, we 
obtained the minimum-energy path in the configuration 
space using the NEB method. Figure [4^i is an illustra- 
tion of the energy profile in the transformation between 
the stable exohedral and endohedral adsorption configu- 
rations, showing the activation barrier (E a ), and reaction 



TABLE II: The activation energy barriers E a and the reaction 
energies E r for the four transitions described in Figure [4] 



System 




E a (eV) E r (eV) 


TR1 - 


► CTR 


2.71 2.37 


NAu~ TR2 - 


> CTR 


2.71 2.22 


TR3 - 


> CTR 


2.74 1.15 


BR3 -> 


In-BR3 


1.51 0.94 




Reactant Transition State Product 



FIG. 4: (a) Schematic energy profile along the reaction co- 
ordinates. The dots represent the energies of the replicas in 
the transformation in the NEB method. The activation en- 
ergy barrier E a and the reaction energy E r are defined by the 
encapsulation barrier from the reactant and the energy differ- 
ence between the reactant and the product, respectively, (b) 
Atomistic model structures for various reaction paths from 
react ants to products. The first three figures display the tran- 
sitions of the N atom from three triangular sites (TR1, TR2, 
and TR3) to the CTR site, whereas the last one depicts that 
from BR3 to In-BR3, the off-centered site. The values marked 
in the middle represent the maximal elongations A^au-Au 
in Afrom the equilibrium distances between two neighboring 
Au atoms located near the N atom during the encapsulation 
process. 



energy (E r ) defined by the energy difference between the 
reactant and the product. Because any triangular site 
can be a gate for atomic insertion, we first considered 
three paths: TR1 CTR, TR2 H> CTR, and TR3 
CTR. In these cases, three Au-Au bonds are maximally 
stretched in the transition states during insertion of the 
dopant atom. We found that the activation energy bar- 
riers of the N atom along all the three paths are almost 
the same (« 2.7 eV). These similar activation barriers 
are strongly related to the maximally elongated Au-Au 
bond lengths, which are essentially the same in all three 
cases, as displayed in the middle column of Figure [4]o 
representing the transition states. On the other hand, 
the reaction energy of the TR3 — » CTR path is about 
1.2 eV smaller than the other two paths, since TR3 has a 
smaller binding energy with a longer Au-N distance than 
TR1 and TR2 (see Table [j). We also checked the BR3 
— >• In-BR3 path, since BR3 and In-BR3 are also stable 
exohedral and endohedral configurations. Interestingly, 
for the BR3 — >• In-BR3 path, a single Au-Au bond cor- 
responding to BR3 is broken to form Au-N-Au and at 
the same time the N atom tend to draw two Au atoms 
located at two nearest AT2 sites. In this transition, the 
activation energy barrier is much smaller to be 1.51 eV 
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than those of three paths through the triangular sites. In 
addition, much more stable In-BR3 configuration than 
CTR for the N encapsulation shows much lower reaction 
energy, which is about 0.94 eV. Consequently, we con- 
clude that the transition path occurring through a single 
Au-Au bond stretch may be preferred in the N atom en- 
capsulation process. The calculated activation barriers 
(Eq) and reaction energies (E r ) are summarized in Ta- 
bled 

The encapsulation events through any triangular site 
would not take place even near the melting temperature 
of the golden fullerene, while the insertion event though 
the BR3 site may occur at relatively low temperatures 
without any external aids. According to our estimation, 
the probabilities that such transitions happen are about 
- KT 11 s- 1 for TR1 CTR, and - 1(T 3 s" 1 for BR3 
—> In-BR3 at 600 K. The insertion events would be even 
enhanced with the help of light or electron-beam irra- 
diation as shown in earlier experiments. It was shown 
that electron-beam irradiation has been used to fuse Ceo 
fullerenes inside host carbon nanotubes 30 and the scan- 
ning tunneling microscope has induced a cis-trans confor- 
mation change in the azobenzene molecule [31 . 

Finally, we calculated the IR active vibrational spectra 
of two most stable exohedral doping configurations (TR1 
and BR3), and an endohedral configuration (CRT). We 
believe that different IR spectra could be used to distin- 
guish structures with different adsorption sites. Figure [5] 
shows our computationally simulated IR intensity data 
for the three configurations. For comparison, we also cal- 
culated the IR spectrum for the bare Au^ 6 cage, displayed 
in solid black in each figure. Its strong IR intensities are 
observed in the frequency range below 200 cm -1 , since 
Au is one of heavy elements, whose frequency is inversely 
proportional to y/m, where m is the atomic mass. We 
found that three modes with dominant intensities arise 
at frequencies of 45, 56, and 96 cm -1 corresponding to 
an Au-Au bond stretching mode and several weak modes. 
In the case of the TR1 configuration, the strongest peak 
appears at much higher frequency of ~ 520 cm -1 as well 
as two weak peaks at « 280 cm -1 and ~ 420 cm -1 due 
to the N adsorbate, as shown in Figure |5^i. For the BR3 
configuration, the strongest peak arises at even higher 
frequency of ~ 580 cm -1 with a very weak peak near 
« 520 cm -1 . It means that BR3 is spectroscopically dis- 
tinguishable from TR1. The direction and the relative 
intensity of each peak are denoted by the green arrows. 

On the other hand, the insertion of the N atom at the 
center combines two lowest modes (45 and 56 cm -1 ) of 
the bare golden cage into one mode while keeping the 
mode around 96 cm -1 , as displayed in Figure^. Three 
other peaks are observed near 90, 140 and 200 cm -1 . As 
denoted by the green arrows in Figure^, all of them are 
due to the vibration motion of the N atom encapsulated 
at the center of the golden fullerene. Especially, the peak 
at 200 cm -1 is composed of three IR- active modes. These 
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FIG. 5: The IR-active modes for (a) the most stable exohe- 
dral adsorption structure (TR1), (b) the next most stable ex- 
ohedral adsorption structure (BR3), and (c) an endohedrally 
doped N@AuJ~ 6 (CRT). The IR-active modes for the bare AuJ~ 6 
cage are displayed with solid black in each figure. Detailed IR 
spectra at wavenumbers below 200 cm -1 are shown in the in- 
sets for the exohedral cases. The model structures with green 
arrows are shown at the peaks to denote the directions and 
the intensities of the IR-active modes contained in each peak. 
In (b), the symbol indicates that the nitrogen atom moves 
towards the center of the Au-Au bond. 

highest IR-active modes occur at much lower frequency 
than those of either TR1 or BR3. This is attributed to 
the flat potential surface of the N atom near the center 
of Au^" 6 resulting in a small force constant. 

CONCLUSION 

In summary, we studied the binding effects of a nitro- 
gen atom on the structural and electronic properties of 
an anionic golden fullerene Auj~ 6 using density functional 
theory. The general trend is that the triangular (hollow) 
and bridge sites are preferred for exohedral nitrogen ad- 
sorption. When the N atom is adsorbed on the cage, elec- 
trons are transferred from Au^~ 6 to the N atom, and the 
binding character tends to be ionic. We also searched for 
the minimum-energy transition paths between the sta- 
ble exohedral and endohedral adsorption configurations, 
and obtained the activation energy and reaction energy 
for each path. It was found that the transitions from 
the triangular sites to the center of Auj~ 6 would not take 
place with an activation energy of about 2.7 eV and a 
reaction energy of 1.2 — 2.4 eV implying no such tran- 
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sition without an external aid. On the other hand, we 
found that the transition through a Au-Au bond may 
take place even at ambient temperatures, since the tran- 
sition energy barrier is much lower. The IR spectra for 
various adsorption configurations were also computed. In 
the exohedral doping cases, the dominant IR peaks oc- 
cur at higher frequencies with larger intensities than in 
the endohedral doping cases. Our study may help de- 
sign new types of small gold nanocluster derivatives. On 
the basis of our present work, further studies will be car- 
ried out on the binding properties of nitrogen-containing 
molecules such as N 2 , NO, NO2 and NH 3 on the golden 
cage. 
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